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ABSTRACT 



A large number of double parton scatterings have been recently 
measured by CDF. The double parton scattering process measures 
CTg//, a non perturbative quantity related to the hadronic transverse 
size and with the dimensions of a cross section. The actual value 
measured by CDF is considerably smaller as compared with the 
naive expectation, namely the value of the inelastic non diffractive 
cross section. The small value of Ue/ / may be an effect of the hadron 
structure in transverse plane. We discuss the problem by taking 
into account, at all orders, the two-body parton correlations in the 
many-body parton distributions. 



1. Introduction 

CDF has recently shown evidence of a large number of double parton collisions [1], 
namely events where, in the same inelastic interaction, two different pairs of partons 
scatter independently with large momentum transfer. The inclusive cross section for 
a double parton scattering is written as [2]: 



CD = D2{xA,x'y^;h)a{xA,XB)c{x'j^,x'^)D2{xB,x'^]h)dhdxAdxBdx'j^dx'^ (1.1) 
Jp- 

^{xatXb) is the parton-parton cross section integrated with the cut off p'j:, which is 
the lower threshold to observe final state partons as minijets, x is the momentum 
fraction, A and B are labels to identify the two interacting hadrons. is a function 
of the product &{xa,xb)o'{x'j^,x'^). Actually the two different partonic interactions 
are localized in two regions in transverse space with a size of order (l/p^)^ and at a 
relative distance of the order of the hadronic radius R, in such a way that the two 
partonic interactions add incoherently in the double scattering cross section. The non 
perturbative input in Eq.(l.l) is the two-body parton distribution D2{x,x';h), which 
depends on the fractional momenta of the two partons taking part to the interaction 
and on their relative transverse distance b. The transverse distance b has to be the 
same for the two partons of hadron A and the two partons of hadron B, in order to 
have the alinement which is needed for a double collision to occur. D2 is a dimensional 
quantity and therefore the process introduces a non perturbative scale factor which is 
related to the hadronic transverse size. 

The simplest possibility to consider is the one where the dependence of D2 on the 
different variables is factorized: 

D2{x,x';h) = f,ff{x)f,ff{x')F{h) (1.2) 

/e// is the effective parton distribution, namely the gluon plus 4/9 of the quark and 
anti-quark distributions and F{h) is normalized to one. Multiparton distribution are 
then uncorrelated and D2 does not contain further informations with respect to the 
one-body parton distribution (actually feff) apart form the dependence on b, whose 
origin is the dimensionality of D2 and which gives rise to the scale factor CTeff- In fact 
in this case one may write 

an-^ (1.3) 



with 



J— = [ F^{h)d% (1.4) 
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and 



I feff{xA)feff{xB)^{xA,XB)dXAdXB, 

Pt 



(1.5) 



the single scattering expression of the perturbative QCD parton model. 

Eq.(1.2) is the basic hypothesis underlying the signature of a double parton col- 
lision which one has been looking for in the experimental search[l,3]. The expected 
characteristic feature of a double collision is in fact that it should produce a final state 
analogous to the one obtained by super-posing two single scattering processes. CDF 
measures: 



By looking at the dependence of cJe// on x CDF has been able to verify the cor- 
rectness of the factorization hypothesis in Eq.(1.2). The range of values of x available 
is limited to a; < .2, for the interaction producing a pair of minijets, and to a; < .4 
for the interaction giving rise to a minijet and a photon. In the limited range of val- 
ues of x available, the factorization hypothesis has shown to be consistent with the 
experimental evidence. 

Since the uncorrelation hypothesis, as expressed in Eq.(1.2), is not inconsistent 
with experiment, one can work out the case where all multiparton distributions are 
uncorrelated and one may look for the sum of all multiparton interactions to the 
hadronic inelastic cross section. The subset where all multiple parton collisions are 
disconnected can be easily summed up in the uncorrelated case [4]. The result is the 
semi- hard hadronic cross section an-, which represents the contribution to the hadronic 
inelastic cross section from events with at least one semi-hard partonic interaction. The 
actual expression is 



The integration on the impact parameter of the hadronic collision, /3, gives the di- 
mensionality to the cross section. The argument of the integral has the meaning of a 
Poissonian distribution of multiple semi-hard partonic interactions with average num- 
ber depending on the impact parameter. 

CTeff is related to cth through Eq.(1.6). Indeed the precise relation requires the 
knowledge of the functional dependence of F{P). Still the expression in Eq.(1.6) shows 
that an does not depend on many details. The only relevant quantity is the value of 
P which makes F{P) small, in such a way that the overall argument of the integral 
becomes close to zero. In the simplest case where F{P) = exp(— /?^/i?^)/7ri?^ one 
obtains a closed analytic expression for an- Actually 



an = 




(1.6) 



an = 27ri?^ [7 + Inn + Ei{k 



(1.7) 
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where 7 = 0.5772 ... is Euler's constant, k = 0-5/(277 i?^) and Ei{x) is the exponential 
integral. For small k one obtains an 2t{R'^k = as, for large k, namely as cxd, one 
obtains an ^ 27ri?^ (7 + /n^) . In this simplest example o"e// = ^tvR"^. The value of an 
is therefore proportional to the measured value of agff, the proportionality factor is 
slightly dependent on energy and on the cutoff. Sensible values of the hadron-hadron 
cm. energy and of the cutoff give values for an which are some 30 — 40% larger with 
respect to the value of aeff. Although an has no reason in principle to be close to ainei, 
a value of an as small as aineif^, which would result from such considerations, seems 
nevertheless too small. Indeed an, as given in Eq.(1.6), has a smooth dependence on 
the infrared cutoff. If one could really have the interaction under control also at small 
values of one would say that an has to be the same as the inelastic non diffractive 
cross section ainei- It would therefore be rather natural to expect for an a value not 
much smaller as compared to the value of ainei- The large difference between an and 
(^inei could therefore be an indication that the simplest assumptions underlying the 
derivation of the expression in Eq.(1.6) have to be revised. 

The main hypothesis which has been done to obtain the expression for an in 
Eq.(1.6) is the Poissonian multiparton distribution. On the other hand one has to 
expect correlations between partons as a consequence of the binding force. While most 
probably correlations will affect the x dependence of the multiparton distribution only 
for finite values of x, and therefore at large rapidities, correlations in the transverse 
parton coordinates are present in every kinematical regime. Indeed the main reason of 
interest in multiple parton collisions, besides the identification of the process itself, is 
precisely the measure of the many-body parton correlations, which is an information 
on the hadron structure independent on the one-body parton distributions usually 
considered in hard processes. 

In the present paper we work out the most general expression for the semihard 
cross section an, which one obtains by assuming that only two-body parton correla- 
tions are present in the many-body parton distributions and by summing all discon- 
nected multiple parton interactions. The two body parton correlation in transverse 
plane is then worked out in great detail in a few cases by considering different explicit 
shapes. 



2. General formalism and semi-hard cross section 

At a given resolution, provided by the cut off p^*"^ that defines the lower threshold 
for the production of minijets, one can find the hadron in various partonic configura- 
tions. The probability of an exclusive n-parton distribution, namely the probability 
to find the hadron in a configuration with n-partons, is denoted by Wn{ui . . . 
Ui = {hi,Xi) represents the transverse partonic coordinate and longitudinal frac- 
tional momentum Xi while color and flavour variables are not considered explicitly. The 
distributions are symmetric in the variables Ui. One defines the generating functional 
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of the multiparton distributions as: 



-^['^]=X^~T / J{ui) ■ ■ ■J{Un)Wn{ui . . .Un)dui . . .dUn, (2.1) 

where the dependence on the infrared cutoff p"^^"^ is implicitly understood, and one 
may introduce also the logarithm of the generating functional: J-'[J] — ln(.Z[J]). The 
conservation of the probability yields the overall normalization condition 



Z[l] = 1. 



(2.2) 



One may use the generating functional to derive the many body densities, i.e. the 
inclusive distributions Dn{ui . . . Un)' 



Di{u) = Wi{u) + I W2{u,u')du' + W3{u,u',u")du'du" + ... 



SZ 
5J{u) 



D2{ui,U2) = W2{ui, U2) + j Wz{ui, U2, u')du' + ^ y Wa{ui,U2, u' , u")du'du" . 



d^Z 



6J{ui)5J{u2) 



52^ 



j=i SJ{ui)5J{u2) 



J=l 5J(ui)5J(u2) 



J=l 



(2.3) 

The many body parton correlations are defined by expanding J-\J\ in the vicinity of 
J= 1: 



T\J\= J D{u)[J{u)-l]du+f^^J Cn{u^...Un)[J{ui)-l] 

n=2 



(2.4). 



[j{un) — 1] dui . . . dUj 



Here D = Di and the correlations Cn describe how much the distribution deviates 
from a Poisson distribution, which corresponds in fact to = 0, n > 2. 

By assuming the validity of the AGK cutting rules [5] for semi-hard interactions 
one may express the semi-hard cross section as the sum of all probabilities of multiple 
parton collisions: 



^H = J d^P<JH{P) 



with 



4 



n 



El S S , f, 

^ ^. SJ'{u[ - (3) ' " 6J'{u'^ - P) ""^^ (2-5) 

n m 

X {i - n n [1 - } n 



J=J'=0 



where /? is the impact parameter between the two interacting hadrons A and B and &ij 
is the elementary probabihty for parton i (of ^) to have a hard interaction with parton 
j (of B) . The semi- hard cross section is constructed summing over all possible partonic 
configurations of the two interacting hadrons (the sums over n and m) and, for each 
configuration with n partons from A and m partons from B, summing over all possible 
multiple partonic interactions. This last sum is constructed asking for the probability 
of no interaction between the two configurations ( actually nr=inj=i[-'- " )■ 
One minus the probability of no interaction is equal to the sum over all semi-hard 
interaction probabilities. 

The presence of multiple parton interactions is induced by the large flux of partons 
which is effective at large energies. The most important contribution to the semi- 
hard cross section, as a consequence, is the contribution of the disconnected partonic 
collisions, namely the interactions where each parton undergoes at most one semi- 
hard collision. These are, in fact, those multiple partonic interactions that, at a given 
number of partonic collisions, maximize the parton flux. We simplify therefore the 
problem by expanding the interaction probability ( the factor in curly brackets ) as 
sums and by removing all the addenda containing repeated indices: 

n,m ^ 

|l_ JJ[1_^..]|^ E ^ij^kl + ... (2.6) 

i,j ij ' ij kyti,lytj 

as a result the semi-hard cross section is constructed with multiple disconnected parton 
collisions only, where disconnected refers to the perturbative component of the inter- 
action. Some features of semi-hard parton rescatterings are presented in Appendix A. 
Because of the symmetry of the derivative operators in Eq.(2.5) one can replace the 
expression in Eq.(2.6) with: 

nmdii — ^^('^ ~ l)m(m — l)a"iia"22 + • • • 

in such a way that the sums over m and n can be performed explicitly. As a consequence 
the cross section at fixed impact parameter, (Th{P), can be expressed by the operatorial 
form: 

(Th{P)= [l - exp{-6 ■ a ■ 6')'^Za[J + 1]Zb[J' + 1] (2-7) 



5 



We have avoided writing explicitly the variables u and u' and the functional derivative 
5/5J{ui) has been simply indicated as 5i. 

The form of aniP) given by Eq.(2.6) is still too complicated to be worked out in its 
general form, since all possible multi-parton correlations are present in Z. Therefore 
we further simplify the problem by taking into account two-body parton correlations 
only. Our explicit expression for is therefore: 

1] = J DA,B{u)J{u)du+^ j CAA'^^^)J{^)J{^)d^d^ (2-8) 

where D{u) is the average number of partons and C(u, f) is the two-body parton 
correlation. 

Either by using techniques of functional integration or by means of a suitable 
diagrammatic expansion [6] one is able to obtain in this case a closed expression for 

aniP) = 1 - exp [-1 En «n - I En bn/n] (2.9) 
where an and bn are functions of the impact parameter P and are given by 



n= J DA{ui)d-{ui,u[)CB{u[ - P, u'2 - P)d-{U2,U2)CA{U2, U3) . . . 

. . . Db (itn ~ /^) H ^'^i^'^'i 



(2.10) 



I CA{un,ui)a{ui,u[)CB{u[ - P,U2 - P)a{u2,U2) . . . 

J (2.11) 

...Cb«_i - P, u'^ - p)a(u'n,Un)Y[duidu[ . 

The actual expression for a„ (that will be referred as "open chain") holds for n 
odd. When n is odd one may also have the symmetric case, where the expression 
begins with Db and ends with Da- When n is even the initial and final distribution 
are either both Da or both Db- It^ the definition of bn, that will be referred as "closed 
chain" n is always even, so that one of the ends is A and the other is B. Sometimes 
the expression n = 2m will be used. One may notice that, at a given order in the 
number of partonic interactions, one can obtain a term of kind a from a term of kind 
b by replacing one C with a pair of D's. The operation can be done in n ways. The 
combinatorial meaning of the 1/n factor multiplying each term of kind b in Eq.(2.9) 
is then understood. The factor 1/2 in Eq.(2.9) is the consequence of the symmetry 
between A and B. 

The cross section is given by an integral on the impact parameter of the interaction 
probability, aniP), that is expressed as one minus the probability of no interaction. 
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The probability of no interaction is given by the negative exponential of the sum over 
all possible different connected structures, namely all structures of kind a„ and of 
kind bn- With our approximations, Eq.(2.6) and Eq.(2.8), these are in fact all possible 
connected structures which can be built with the average numbers Da,b, the two-body 
correlations Ca,b and the interaction a . Expanding the exponential, the cross section 
can then be expressed as the sum over all possible structures, both connected and 
disconnected. 

One will notice that, when no correlations are present, all terms of kind b disappear 
and only the first of the terms of kind a, namely Dao^Db is left. In that limit the 
cross section is given simply by: 



d^P<l- exp 



Da{u — P)DBiu')a{u,u')dudu' 



(2.12) 



which corresponds to the Poissonian distribution discussed in the introduction. 



3. An explicit case: Gaussian correlation 

The bulk of the semi-hard cross section originates in the region of small x values. 
The experimental observation on the independence of o"e// as a function of x, at small 
x-values, on the other hand indicates that correlations in x, at small x, are not a strong 
effect in the two-body parton distribution. A sensible case to consider is therefore the 
one where there are no correlations in fractional momenta, nor between fractional 
momenta and transverse coordinates, while only transverse coordinates are correlated: 

Diu) = fix)Dih) , Ciu, u') = fix)fix')Cih, W). (3.1) 

Given the localization of the partonic collisions in transverse space, the dependence of 
the parton-parton interaction probability on b and b' is represented as a 5-function: 

a{u,u') = aa,,^>S{h-W). (3.2) 

All integrations on the fractional momenta x and x' are then factorized from the 
integrations on the transverse coordinates and result in the single scattering term cross 
section as, whose expression is given in Eq.(1.5). an and bn are therefore considerably 
simplified: 

an = a^ [ DA{h^)CB{hi-p,h2-P)CAih2M)---DBihn-P)l[dhi (3.3) 



bn = (T^ J CA{hnM)CB{hi-pM-P)---CB{hn-l-p,hn-p)lldhi.. (3.4) 
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To proceed further one needs to consider explicit functional forms for C(b, b'). A 
simplest extension of the Gaussian model discussed in the introduction corresponds to 
the following Gaussian expression of the correlation term: 



L»(b) = /iexp[-Fb2] , C(b,b') = A;exp[-|F(b + b')^-^G'(b-b')^]. (3.5) 
After defining b = y + 6^ is explicitly written as: 



bn = a^sk^'J exp[-|nF/32] exp[-(F + G) Ey.'] exp[-(F -G)Zy^■ y^+l]Udy^, 



with the convention that the variable yn+i coincides with yi . 
bn can then be worked out through the substitution 



which leads to 



(3.6) 



bn = Ja^k^ exp[-inF/32] / exp[- ^ v,'] 11 

where the Jacobian, as discussed in Appendix B, is J = (//" — . bn is therefore 

finally written as: 



bn = a^kViijJ' - zy")-2 exp[-inF/32 



(3.60 



The case of the open chain, an-, is less symmetric and requires a slightly less 
straightforward treatment. Using the y variables, defined previously, the expression 
of an takes the form: 



= a^/i^/c-i exp[-inF/32] /exp[-(F + G)y^ -M- y]]\dyi, (3.7) 



where 



and 



M = 



y'^ = (yi,y2,...,yn) 



/l + r r ... 
r 1 r ... 



(3.8) 



\ ... r l + r/ 

r = i(F-G)/(F + G). 



(3.9) 



(3.10) 
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The Gaussian integral can be evaluated (details of the calculations are reported in 
Appendix B) leading to the relation: 



an = a^h^k''-^Tv''r''-^[{2r + l)Un-i{l/2r)]-^ exp[-^nF(3^] ; (3.7') 

where Un is the Chebyshev polynomial of second kind [7]. 

While the terms and bn are computed exactly, the sum of the series in Eq.(2.9) 
can be performed only in limiting cases. We discuss the case where the correlation 
length is much smaller with respect to the hadronic radius. The parameter G char- 
acterizes the correlation while F is related to the hadronic transverse size. Small 
correlation lengths correspond to F << G. We work out therefore the leading order 
term in F while keeping the full structure in G both in and in In order to 
find the limiting expression of when F « G, one needs to work out the limit of 
Un{l/2r) for l/2r ^ -1 + e. One obtains: 

C/„(-l + e)^«(-)-(n+l) 

By using this expression for Un and with similar, but simpler, manipulations one 
obtains the following limiting form for a„: 

(3.7") 



4:nkF 



with 

^)a5exp[-iF/?^] (3.11) 
The whole series of Eq.(3.1) can be then be summed yielding as a result: 

Sa=^[-Hl-Z)]. (3.12) 

To work out the limiting case for bn we keep the exact value of the difference 
{/J, — v)'^ = 2F, while setting ji ^ v ^ ■\/G/2 everywhere else. Keeping moreover into 
account the condition n = 2m, we obtain: 

The resulting sum of the series in Eq.(3.1) is therefore given by 

St = ^MZ') ■ (3.13) 

where Cs{x) = 'Y^x^ /n^ is the trilogarithm function[8]. 

The final limiting form for the cross section, when F << G, is therefore: 
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an = / d'^^aHiP) = / d''p{l - expf-JS'a - \Sb]} , (3.14) 



In the limit G oo the correlation goes to zero. In this case one obtains Z ^ 0, 
Sa ~^ {nasf^^kF) exp[— ^F/?^] and — > 0. The expression of an in Eq.(3.14) is then 
reduced to the Poissonian distribution discussed in the introduction. 

4. Two qualitatively different features of the correlation term 

The uncorrelated multiparton distributions are characterized by two independent 
features. One feature is that the n-body parton distribution is factorized in the product 
of n times the one-body parton distribution. The second feature is that the distribution 
in the number of partons (namely after integrating over all other degrees of freedom) 
is a Poissonian. Both features are affected by the introduction of the two-body corre- 
lations and it is therefore interesting to study the two effects separately. 

One may modify the number distribution, without introducing non-factorized 
two-body correlations in b, by using the factorized expression 

C{u,u') = ^\D{u)D{u') (4.1) 

The terms appearing in Eq.s (2.10,11) are very much simplified in this case. After 
integrating over the longitudinal variables, the transverse integrations always appears 
in the unique form 

T = T{I3) = asj DA{h)DB{h - P)dh (4.2) 

One therefore obtains: 

|En«n=T/(l + AT) 

^E^&2n^/2m = ln(l-A2T2). 
and the cross section, at a fixed value of the impact parameter, is expressed as: 

aH{(3) = 1 - [1 - A2T2]-V2exp[-T/(1 + AT)] . (4.3) 

The parameter A represents a measure of the deviation from the Poissonian distribu- 
tion, both for the initial states and for the elementary collisions. In Appendix C some 
further elaboration on the factorized case are presented. 

Perhaps more interesting is the case where the number distribution of the in- 
coming parton is Poissonian, because the correlation integrates to zero, while the 
resulting distribution in the number of parton collisions at fixed impact parameter is 
non-Poissonian, as a consequence of the presence of the correlation term. To study 
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this case a different way of dealing with the approximation introduced in obtaining 
Eq.s (3.6", 7") is useful. 

The alternative procedure is the following: The variables are substituted by 
the mean value Y = ^y/n and by the differences Xj = yi+i —yi, with the constraint 
^ Xj = 0. The relevant Jacobian is J = 1 and so the integration volume is transformed 
according to: H'^Yi = cfV (ixj(5(5^ Xj). Then in the terms containing F, which 
simply defines the size of the hadron, one performs the substitutions yi ^Y/n whereas 
the terms containing G and the differences x^ are not modified. For the calculation 
of ttn, since only n — 1 differences appear, the integrations in dY and c/xj decouple 
in the Gaussian approximation and the result of Eq.s (3. 7", 11) is recovered; for the 
calculation of bn the constraint on the x^ is essential, it may be implemented in the 
standard exponential form S{J2^i) = /(Stt)"^ exp(i ^ q • xjdq and carrying out all 
the integrations the result of Eq. (3.6") is reproduced. 

This alternative calculation suggests an initial distribution which can be easily 
worked out and which yields to a pure Poissonian distribution in the number of partons 
in the initial states. One may in fact consider the following expressions: 



D{h) ^ hexp[-Fh^] , C(b,b') = -/c'A2exp[-iF(b + b')2 - ^^(b-b')^] (4.4) 



where A2 indicates the two-dimensional Laplace operator acting on the transverse 
difference b — b' or y — y'. By integrating C over the transverse variables one obtains 
zero. The number distribution depends therefore only on D and it is Poissonian by 
construction. The parameter k' is obviously related to the parameter k in the previous 
paragraph, k was the strength of the correlation at b = b' = 0. To keep the same 
normalization we make the position k' = k/G. 

The contribution of all open chains with n > 1 is therefore zero. As previously 
observed in this case the differences have no constraint. Every independent integra- 
tion over a transverse difference corresponds therefore to an integration of a single 
correlation term, which gives zero. Small nonzero effects have been possibly lost in 
the approximate substitution Y/n in the terms containing F in Eq.'s(3.6',7'). 

When one computes the contribution of the closed chain bn, the constraint is ef- 
fective. The calculation is again performed by means of the exponential representation 
of the 5-function and the resulting expression can be brought to the following form 



The final result is obtained by summing all the terms above, with the weight (l/2m). 
The cross section (at P fixed) is therefore 



G 



(Z/m)2"^(2m)! 



(4.5) 





(4.6) 



m=l 
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While a closed analytic form for the sum Sf, is not easily written in the general 
case, in the limiting case where the number of partonic collisions is large one may use 
the Stirling approximation for the factorial: m~"^ ~ V27rme~"^/m!. With the help of 
the decomposition (2m)! = 2^"^m!(|)^ the sum can be expressed as 

^6 = ^^{w) , (4.7) 

where 



\ 2Jm „„m 

ml 



iIj{w) can be evaluated by computing, as intermediate step. 



' ml 



1 

and by performing the integrations afterwards. The resulting expression is: 

V'(^x;) = 2£2(i[l- v^T^]) -ln'(i[l + ^rT^) . (4.8) 
where C2{x) — jn^ is the dilogarithm function. 

5. Concluding discussion 

The small value of <yeff^ the dimensional parameter characterizing double parton 
scatterings, which has been measured recently by CDF, is an indication that two-body 
parton correlations, in the many-body parton distribution of the proton, are likely 
to be sizeable. In the case of an uncorrelated many-body parton distribution, the 
value of (Jeff puts a constraint on the range of possible values of a^/, the semi-hard 
contribution to the hadronic inelastic cross section. The actual measured value of (Jeff 
would give rise to values of ou of the order of Oi^eil'^ also at very large cm. energies, 
where one would rather expect (Th — (Jinel- The experimental evidence is also that, 
in the x region accessible experimentally namely at small x values, the correlation in 
fractional momenta is not a large effect. 

In the present paper we have worked out the semi-hard cross section an in the 
case where initial state partons are correlated in the transverse parton coordinates. In 
fact an can be worked out rather explicitly when only two-body parton correlations 
are included in the many-body parton distributions and when each parton can have at 
most one semi-hard interaction. There are two qualitatively different features in the 
two-body parton correlation, and both change the relation between an and a"e// with 
respect to the uncorrelated case: 
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1- The distribution in the number of partons is not any more a Poissonian, although 
the dependence on the kinematical variables of the different partons is factorized. 

2- The overall distribution in the number of partons, namely after integrating on 
the partonic kinematical variables, is a Poissonian but the dependence on the 
partonic kinematical variables is not factorized, in which case the two-body parton 
correlation integrates to zero. 

The general case is obviously a combination of the two possibilities. We point out 
however that both cases separately can give rise to a small value of a^f / while keeping 
the value of an close to (Tinei- 

In the first case an is obtained by integrating the expression in Eq.(4.3). In the 
second by integrating the expression in Eq.(4.6). The critical value of the impact 
parameter /3c, which gives the size to the cross section ajj, is the value which makes 
small the argument of the exponential in the expression of aniP)- The detailed de- 
pendence of the argument of the exponential ai (3 < (3c is not of great importance for 
the determination of an- It only determines the degree of blackness, with respect to 
semi-hard interactions, of the overlapping hadronic matter at the given value of the 
impact parameter. It is therefore enough to determine the argument of the exponen- 
tial at the edge of the interaction region. In case 1 the behaviour of the argument of 
the exponential when going at the edge of the interaction region corresponds to the 
limiting behaviour of the expression below at small T values: 

+ lln(l - A^T^) ^T = asJ DA{h)DB{h - P)d% 

while in case 2 it corresponds to neglect Si, as compared with T in Eq.(4.6): 

T + Sb^T = as J DA{h)DB{h - (3)d%. 

In both cases the critical value f3c has the same value which one finds for the uncor- 
related distribution and therefore uh, as in the uncorrelated case, is roughly equal to 
27r/F, where F is the parameter giving the extension of the one-body parton distri- 
bution in the transverse plane. 

When the double scattering cross section is worked out by expanding the expres- 
sion of one finds for the effective cross section the relation 

27r 

""'^^ ~~ F-(l + A)2 

in case 1, and 

= [ C(b, b')C(b -P,h' - P)d^bd^b'(fp 

aeff 27r J 

in case 2. 
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A qualitative feature is that in both cases one obtains a value of erg// which may 
be sizably smaller with respect to 2ti / F ^ cth- While, on the other hand, nothing 
prevents the value of ch to be close to the value of Cinei- The smaller value of Ueff, 
with respect to the expectation of the uncorrelated case, is rather generally associated 
with the increased dispersion of the distribution in the number of partonic collisions: 
In the case of no correlations the distribution is strictly Poissonian when the impact 
parameter is fixed. When correlations are introduced the distribution in the number of 
parton collisions, at fixed /3, is not Poissonian any more and the natural consequence 
is that the dispersion in the number of collisions is increased. 

The indication from the measure of the rate of double parton scatterings is there- 
fore that two-body parton correlations arc likely to be important while, unfortunately, 
one cannot say much about dynamical quantities, like the the correlation length. Use- 
ful observables to be measured, in order to get some more insight into the problem, 
would be the semi-hard cross section an and the triple parton scattering cross section. 
Our present analysis shows that an can be reliably discussed in perturbation theory. 
The measure of an, in association with (Te//, would help considerably in clarifying the 
size of the effect induced by the presence of the two-body parton correlations: All the 
considerations of the present paper are based on the prejudice that au should have a 
value rather close to the value of ainei- 

The measure of triple parton scattering would allow one to test the possibility discussed 
in case 1, since in that case the rate of triple scatterings would be strictly fixed by 
the measured rate of double scatterings. In the other cases presently discussed the 
knowledge of the rate of triple parton collisions would allow one to obtain the actual 
values of the parameters of the correlation term. 

A lot of effort has been put on the study of the proton structure as a function 
of the momentum fraction x. The distribution of partons however depends on three 
degrees of freedom, the momentum fraction x and the transverse parton coordinate b. 
The measure of the rate of triple and of higher order partonic collisions is the essential 
tool to learn on the parton structure of the proton in transverse plane. 
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Appendix A 



In this appendix some short observations about the effects of the rescattering 
processes wiU be displayed. The general form of the hadron-hadron scattering cross 
section Eq. (2.5) has been reduced to a more manageable form by throwing away all 
the rescattering process by means of the position: 



n,m ^ 



(Al) 



It is possible to deal with some rescattering process without too much effort if we look 
at processes where only one of the two colliding partons has already suffered some 
collision. In this case instead of Eq. (A.l) we write 

n,Tn ^ 



1-J 



ij ki^i, 



{A.2) 

In this expression the first term shows that the partons belonging to the first hadron 
suffer only one collision, the partons belonging to the second hadron undergo any 
number of collisions, the second term describes the symmetrical situation, the third 
term eliminates the double counting of the single-collision processes. *. 

The operatorial form of the cross section at fixed impact parameter aniP) acquires 
now the more complicated shape: 



aniP) = [l-exp((5 • [e* '^' - 1]) - exp(5' • [e"^'^ - 1]) + exp(-(5 ■ a ■ 5') 
Za[J + 1]Zb[J' + 1 



X 



(A.3) 



J=J'=0 



This expression seems of uneasy interpretation but is gives some information when re- 
duced to simpler particular cases The simplest, but non trivial result is produced when 
the incoming parton distribution is purely Poissonian and only the double scattering 
is actually considered. In these situations it is easily seen that instead of Eq. (2.10) 
one obtains the form: 



=/*{ 



1 — exp 



Da {u)a{u, u')DB{u')dudu' 



X 



exp 



I / DA{u)a{u,u')a{u,u")DB{u')DB{u")dudu'du' 
\ J DA{u)DA{u")a{u,u')a{u",u')DB{u')dudu'du" 



+ 



- 1 



{AA) 



This possibility was already considered in ref [6], for an unsymmetrical situation 
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This form is not unexpected, it could have been written by hand, Eq.(A.3) shows a 
more systematic way of deriving it and possible further corrections. The final distri- 
bution of the collision is not Poissonian and this deviation is purely due to the hard 
dynamics. It can be observed that precisely this origin in the hard collision may also 
offer a way to distinguish these eflFects from the effects of correlations in the incoming 
two-body distributions, in fact in the case here sketched there should be an unbal- 
ance among the kinematical variables of the pair of jets originating from the scattered 
partons. 



Appendix B 



The Jacobian appearing in Eq.(3.6), which arises from the transformation from 
the variables y to the variables v, is expressed as J = [detj']~^. The actual form of 
J is: 

/ IJ -u ... 0\ 
n -ly ... 
J= . . : : 

\-u ... fi/ 

and the result detj^ = ii^ — u^ is easily obtained by expansion according to the last row. 
The exponent —2 appears because the matrix describes the inverse transformation and 
there are two transverse dimensions. 

In order to calculate the determinant of the matrix M. in Eq.(3.7) an auxiliary 
matrix A is introduced so that we have 



M 



/I + r r 
r 1 r 

V 








\ 




1 + r/ 



A 



/I r 

r 1 r 



0\ 




\0 



1/ 



By means of the standard rules for the computation of the determinants two 
recurrence relations are obtained: 

detM-n = detAn + 2r detAn-i + r^detAnz 

detAn = detAn-i — r'^detAnz 
which imply the simpler relation 

detMn = (2r -|- l)detAn-i 

The recurrence relations for the determinants of An are very similar to the recurrence 
relation for the Chebyshev polynomials Un [7] and in this way it is possible to get the 
final expression for them 

detAn = r''Un{l/2r). 
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For completeness we give also the explicit expression of the ?7-polynomials: 

Unicose) = — . 



Appendix C 

In this appendix we shall briefly discuss some non-Poissonian one-body densities, 
the essential point were already shown in ref.[6], but it may be useful to state again 
them in order to have a comparison of the previous treatment. The starting point of 
Eq.(2.7) is specialized to the case of pure one-body densities in the form 

an = (1 - exi^Tdd')XA{y)XB{y)\y=y„,y'=y'„. {C.l) 

The Poissonian distribution is given hy X = ex.p[y — the situation described in 
Eq.(4.1) corresponds to 

X = exp[y - yo] + ^x[y - yo]^ , {c.2) 

Then by applying Eq.(C.l) to Eq.(C.2), one gets as intermediate step 

exp[-Ta' + ^X{Td'r] exp[y - X] + |A[y - y^^ly^y^, {C.3) 
By the use of the general formula: 

1 r qx'^ 



ex.p[-pdj exp[qx ] = exp 

Vl - 4pg 



1 + Apq 



it is seen that the expression appearing in Eq.(C.3) reduces to the form already given 
in Eq.(4.3). 

One could also choose the form 



x = [i-{y-y,)]-^ 

which corresponds to a negative binomial distribution for the initial partons and obtain 
in this way the result 

aH(/?) = l-T-"t/[a;l;l/T], 

by representing the incoming distributions as Laplace transforms as shown in ref.[6], 
note that a slight simpliflcation has been introduced in the distribution, with respect 
to that reference, and this reflects into a slight simpliflcation of the result, which keep 
however its main properties, in particular U is the irregular confluent hypergeometric 
function. 
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